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Abstract. We discuss the computation of the three loop anomalous dimensions for various 
operators used in deep inelastic scattering in the MS and RI' schemes. In particular the results 
for the n = 5 and 6 Wilson operators in arbitrary linear covariant gauge in the RI' scheme are 
new. 
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1 Introduction. 



The use of symbolic manipulation and computer algebra has been an invaluable tool for provid- 
ing large loop order results in quantum field theory in general and in quantum chromo dynamics 
(QCD) in particular. The outstanding example is that of the full three loop anomalous dimen- 
sions for flavour non-singlet and singlet unpolarized operators for deep inelastic scattering as an 
analytic function of the operator moment n, [H El El H] . Also the Wilson coefficients have been 
provided to the same precision. This large project, of the order of ten years, required not only 
the extensive use of the symbolic manipulation programme Form, [5], but also its own devel- 
opment to handle the unforeseen complexity of the computation. An earlier approach, [61 [7], 
to this problem in deep inelastic scattering was to determine results for fixed (even) moments 
using the Mincer algorithm, [8j, translated into FORM, [9]. For instance, the first few moments 
were given in [6], [7j and subsequently those for n = 12 and n = 16 appeared in [101 lllj. At that 
time the fixed moment expressions were used to parametrize the full expressions but such an 
approach was clearly incomplete lacking the correctness of a full evaluation. However, the re- 
sults subsequently served as very important independent checks on the final arbitrary n results. 
Now that the computational algorithm has been established, in principle it can be applied to 
other operators underlying related phenomenology. For example, the case of polarized Wilson 
operators will be relevant for spin physics. In addition in the spin context there is interest 
in a similar operator called transversity, [121 [T3l 114] . This corresponds to the probability of 
finding a quark in a transversely polarized nucleon polarized parallel to the nucleon versus that 
of the nucleon in the antiparallel polarization. From a theoretical point of view it is similar 
to the non-singlet unpolarized Wilson operator but experimentally it is not as accessible since 
there is no direct coupling to quarks. Nevertheless there have been proposals to study it at 
RHIC. Therefore, whilst in principle it is possible to calculate the arbitrary moment three loop 
transversity operator anomalous dimensions in the MS scheme, it would be important to have 
strong independent checks on any future full result. Akin to the 1990's approach for the Wilson 
operator there is therefore a need for fixed moment calculations. Aside from this motivation, 
there is a secondary one. 

One of the ingredients necessary to study the structure functions is the measurement of 
the non-perturbative matrix elements. From the theory point of view, a tool which achieves 
this is lattice regularization and various groups, such as QCDSF, have developed a substantial 
programme to determine key matrix elements. (See, for instance, |15| . I16 (. H7].) However, one 
technical aspect of such work is ensuring that the results agree in the continuum with expec- 
tations from the ultraviolet limit. One approach in this respect is for the lattice results to 
be matched onto the perturbative expressions in the chiral limit, where to aid precision one 
would prefer the results to as high a loop order as is calculationally feasible. This has been 
considered in a series of articles, [TH1 IHl EO] , to three loops. One technical issue is that to keep 
time (and money) to a minimum, the lattice computations are performed in renormalization 
schemes known as regularization invariant (RI) and its modification, RI', |21| I22j. Unlike the 
MS scheme, they are mass dependent renormalization schemes. Results in this scheme have then 
to be converted to the standard MS scheme. In the continuum QCD has been renormalized at 
high loop order in both RI and RI' in [231 [T8] and the conversion functions established for various 
quantities of interest. Therefore to aid lattice computations of the matrix elements one requires 
the finite part of the analogous Green's functions but only for low moment since clear signals for 
higher moments information are hard to extract from the numerical noise on the lattice. Given 
the need for such accurate results for a specific Green's function on the lattice, not only for the 
transversity operators but also for the Wilson operators, we report the results of recent com- 
putations in this area, [201 124] . One consequence of the Green's function considered is that the 
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operator anomalous dimensions emerge as a corollary in the RI' and MS schemes. In addition 
to the results already available, [181 032 1201 121]) we will provide the RI' anomalous dimensions 
for the n = 5 and 6 Wilson operators in arbitrary covariant gauge. Although the finite parts are 
required for lattice calculations for low moment, our computations have been extended to high 
moment for the transversity operator but only to determine the anomalous dimensions. At the 
appropriate point we will mention some of the computer algebra aspects of this and previous 
work which without the power of Form would have rendered the determination of any result 
virtually impossible. 

The paper is organised as follows. In section two we introduce our notation and compu- 
tational strategy in more detail before discussing the appropriate points of the RI' scheme in 
section three. Section four details a simple low moment example, whilst the symbolic manipu- 
lation issues are recorded in the context of the higher moment calculation in section five. The 
explicit anomalous dimensions are given in section six with a few concluding remarks provided 
in the final section. 



2 Background. 

The two basic classes of operators we will consider are the non-singlet Wilson operators 

0^- Un = S%^ i 'f 1 D V2 ...D^ilP (2.1) 

and the transversity operators 

p...u n = S fcj^D v *...D Vn i$ (2.2) 

where ip l is the quark field, 1 < i < Nf for Nf quark flavours, is the covariant derivative 
and a^ u = ^[7^, 7^]- The operation of totally symmetrizing with respect to the Lorentz indices 
{y{\ and ensuring the operator is traceless is denoted by S where the respective but different 
tracelessness conditions are given by 

V ViVj 6w" Vi "- Vi '" Vn = (2.3) 

and, [25], 

^Op-^n = (i > 2 ) , // i , F ,CV / "' -" = . (2.4) 

For (|2,2p the anomalous dimensions are available for all moments n at one and two loops in the 
MS scheme in [261 l27l l28~l [25| [29] . At three loops fixed moment results are available for moments 
up to n = 8 for MS and n = 7 for RI', [191 US 121]. The specific Green's functions relevant for 
the lattice matching are 

G^~ u "(p) = (V(-P) 0%" M »(P) #>)) 
G^- v -(p) = (4>(-p)Op- u "(0)j(p)) (2.5) 

where the operator is inserted at zero momentum. This allows for the application of the Mincer 
algorithm, [3(9], which determines the finite part of scalar massless two point functions using 
dimensional regularization in d = 4 — 2e dimensions to three loops. Unlike earlier approaches 
to extract anomalous dimensions we do not contract the free Lorentz indices with a null vector. 
Whilst such a contraction has the effect of excluding the trace terms in the operator itself or 
the Green's function decomposition, the main reason why we cannot follow that route here is 
that the lattice makes measurements in different directions of the momentum components. This 
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allows for the extraction of the values of each of the individual amplitudes into which the Green's 
functions are decomposed. 

From the form of the operators there will be n n-point Feynman rules for both the Wilson 
and transversity operators, each with two quark legs. However, at the three loop order we will 
work at, only the Feynman rules up to and including three gluon legs will be necessary. Hence, 
there are 3 one loop, 37 two loop and 684 three loop Feynman diagrams contributing to each 
Green's function. These are generated electronically using the Qgraf package, [30], before 
being converted into Form input notation to allow for the application of the Form version of 
the Mincer package, [U [9]. As the Mincer algorithm is only applicable to scalar Feynman 
integrals, for each moment n the Green's functions, (12. 5ft . need to be decomposed into invariant 
amplitudes and Lorentz tensors which respect all symmetry structures. Whilst we will discuss 
more details later, it suffices to note at this point that for Wilson operators there will be two 
independent amplitudes but three for the transversity case. 

Although the lattice computations are ultimately performed in the Landau gauge, we will 
compute in an arbitrary linear covariant gauge. The associated gauge parameter will act as an 
internal checking parameter since, for instance, in the MS scheme the anomalous dimension of 
gauge invariant operators are independent of the gauge parameter, [311 [32]. As the computations 
are clearly automatic we employ the procedure of |33j where all the diagrams are computed for 
bare coupling, g Q , and gauge parameter a . To extract the anomalous dimension the renormal- 
ization constant is fixed (either in the MS or RI' schemes) by rescaling these variables by the 
known coupling constant and gauge parameter renormalization constants 

9o = Z g g , a Q = Z~ x Zaol (2.6) 

in our conventions, where Za is the gluon wave function renormalization constant. The remain- 
ing divergence for each of the Green's function is absorbed into the operator renormalization 
constant together with a specified finite part in the case of the RI' scheme to leave the finite 
parts for the lattice matching. In practice the results are determined in the MS scheme first, 
primarily due to more consistency checks being available before extracting the RI' expressions. 

3 RI' scheme. 

We briefly review the parts of the RI' scheme needed for the computations discussed here. 
Originally the scheme was invented in the context of the lattice, |21|. I22j. but it is not restricted 
to a discrete spacetime. The continuum analogue has been studied to three and four loop 
order in |23|. 118] . In general terms it is a mass dependent renormalization scheme where the 
renormalization of the quark field is chosen to be non-minimal in a way which is appropriate for 
lattice analyses. The coupling constant (and thus vertex) renormalization is performed in an MS 
way and so in some sense the RI set of schemes sits between MS and MOM type schemes. To 
reduce time (and cost), since taking a derivative on the lattice requires significant computation, 
the quark 2-point function, Sw,(p), is renormalized according to the RI' prescription, |2H 122]. 

= i> (3.i) 

p2 — ^2 

where /i 2 is the renormalization point. In other words there are no 0{g 2 ) corrections to E^(p) 
after renormalization at p 2 = /i 2 as these finite parts are absorbed into the quark wave function 
renormalization constant, Z'^. We use the notation throughout that a ' on a quantity indicates 
that renormalization has been performed in the RI' scheme. Otherwise the scheme is MS. For 



lim 

e-+0 
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completeness the RI scheme, which is not of interest to us here, involves taking a momentum 
derivative of E^(p) first before choosing the result to be the tree value at the renormalization 
point, |21l I22j . As an extension of the RI' scheme in the continuum, the gluon and Faddeev- 
Popov 2-point functions are renormalized analogous to (|3,ip . However, as most interest in 
general is in quark 2-point Green's functions, there is no real need to pursue this route, unless 
one was perhaps intending to consider supersymmetric theories. Similar to the lattice we are 
ultimately interested in converting results from RI' to MS and therefore the variables in each 
scheme need to be related. Using the standard conversion definitions 



a 



Z a 



a 



Z 'n. 



(3.2) 



where a = g 2 /(16ir 2 



a 



we have the one loop relations, [23j H8] . 
a + 0(a 5 ) 
1 + ((-9a 2 - 18a - 97) C A + 80T F N f ^j 



a 
36 



a + 0(a 2 



(3.3) 



The explicit expressions to three loops are available in [18]. Though it is worth noting that the 
Landau gauge is preserved in changing between RI' and MS. To illustrate the effect the schemes 
have on the basic anomalous dimensions, we note 
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aC F a + - (a 2 + 8a + 25)C A C F 
4 L v 



6C 2 



' 9a 3 + 45a 2 + 223a + 225 ) C A 



70 (a) = aC F a + 

- 54<7 F - (80a + 72) T F N f ] 
where the group theoretic quantities are defined by 



8C F T F N f 



C F a 2 
36 



Tr T a T 



T F 5 



ab 



rparj-ia 



C F I 



jacd jbcd 



C A S' 



ab 



(3.4) 



(3.5) 



for a colour group with generators T a and structure functions f abc . Clearly the difference in the 
numerical structure in (|3.4p only appears at two and higher loops. 

For the flavour non-singlet operators we are interested in here, we follow a similar route to 
(|3.ip for defining the operator renormalization constant in the RI' scheme. Writing (p) as 
the amplitude in the Lorentz decomposition of (12. 5j) which contains the tree, (T), part of the 
operator, we set 



lim 



T 



(3.6) 



where T is the value of the tree term of amplitude, which may not necessarily be unity given 
the specific (non-unique) way of carrying out the decomposition. 



4 Simple example. 

We now illustrate the preceeding remarks by discussing the case of the n = 2 transversity 
operator in more detail, [19]. First, given the symmetry properties (|2.4p the explicit traceless 
symmetrized operator is 

(d-l) 

+ Jd~i) fr V i>° pX D^ + n ^a vX D x ^) (4.1) 
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in d-dimensions. Inserting (|4.1|) into the Green's function Gj? p {p), it is decomposed into the 
three invariant amplitudes as 



G^ p {p) = T^\p){a^pP + a pp p u 



p 

+ ^\p) [^ u a pX p x + r] pp a uX Px - (d + l) V up v pX Px 

(d-l)(d + 2) uA „ . 
+ ^ >-o pX p v p p p x 

+ (p) (a vX p p p p p x + a px p p p v p x + do pX p v p p p x - rf p <y pX p x p 2 ) (4.2) 

in d-dimensions. It is worth noting that this and other decompositions are not unique since one 
can always take a linear combination of the three (independent) tensor structures consistent with 
the symmetry and traceless properties to form another set of independent amplitudes. However, 
with this choice one can algebraically form a scalar object which is computed via Mincer. For 
instance, [19], 

(d + 2) A 
o-^Pp + VupPv — 2 — (y^xPvPpP 



E ^ = - 8(d-l)(d-2) tr 



rtr 



(°v\PpPpP X + <?p\PpPvP X 

+ da fJ/X p u p p p x - y]v P (yp X p x p 2 ) G^ p (p) . (4.3) 



8(d- l){d-2)p- 



This together with S^(p), i = 2 and 3, are the objects of interest for the lattice matching and 
have been determined to 0(a 3 ), |19j . For this specific example, we note that the construction of 
the tensor basis as well as the amplitude decomposition can easily be carried out by hand. This 
is primarily due to the small number of free Lorentz indices present. Clearly for the extraction 
of the anomalous dimensions and amplitudes for the higher moment operators, such a procedure 
would be unacceptably time consuming by hand. Moreover, it would be prone to elementary 
algebraic errors. 



5 Higher moment issues. 



To extract the anomalous dimensions for the higher moment operators, it is clear that one 
has to proceed with a computer algebra construction to determine the basis for the independent 
amplitudes and hence the projections. We discuss the issues in relation to the n = 8 transversity 
operator as an example, [24j. For this case there are initially seventeen potential tensors into 
which the Green's function (|2.5p can be decomposed. These are built from the relevant vectors 
and tensors of the operator in question, which for transversity are p^, rf" v and a^ v . The only 
constraint being that the Green's function has nine independent indices. Given these seventeen 
tensors then within FORM it is straightforward to construct the seventeen tensors which have 
the correct symmetry, but not traceless, properties. Taking a linear combination of these new 
objects with as yet unrelated coefficients, the relationship between these are fixed by imposing 
the remaining traceless criterion. In practical terms we take successive pairs of free indices and 
contract them. The coefficients of the resulting tensors produce constraints on the seventeen 
initial coefficients which can be solved. Whilst there are more contractions than coefficients 
there is redundancy in the system of linear equations which determine the coefficients. This 
is due to the symmetry of the operator itself. However, there is no unique solution and three 
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coefficients remain unrelated producing three independent amplitudes. (For the Wilson operator 
the corresponding number is two.) For the higher moments, as we are ultimately interested in 
the anomalous dimensions, the specific linear combination one uses is not a major issue. The 
only constraint is to choose that projector of the three which leads to the lowest computation 
time when Mincer is applied. The test for this is to compare the run times for each projector 
to do the full two loop calculation in an arbitrary linear covariant gauge before generating the 
results for the three loop diagrams. 

The other main computer algebra issue is the construction of the Feynman rules for each 
operator. Given that the Green's function has free indices one in principle has to construct 
the full symmetrized and traceless operator before applying the Form routine to generate the 
explicit Feynman rules for the operator. However, given that the Green's function will be 
multiplied by a projector which is traceless, that part of the operator containing rj^ tensors 
will automatically give zero upon contracting with the projector. Therefore there is no need to 
have an operator which is traceless; only an operator which is symmetrized will suffice. This 
will reduce computation time since otherwise with a traceless operator there will be an internal 
intermediate expression swell which will be sorted by Form to produce the equivalent scalar 
expression as ignoring the traceless part. For instance, for the n = 8 transversity operator the 
expression swell would have been substantial. Finally, in relation to the Feynman rules, only the 
part of the operator up to and including two quark and three gluon leg insertions are required 
for the full three loop computation. To illustrate the size of our higher moment calculation the 
Form module containing the operator Feynman rule was 36 Mbytes for n = 7 transversity and 
300 Mbytes for n = 8 transversity, [23]. Indeed the latter calculation could only be performed 
in the Feynman gauge rather than the full linear covariant or Landau gauges. Even then it took 
of the order of 40 days on a dual opteron 64 bit SMP 2GHz machine. Hence, only the MS result 
was determined with the RF scheme anomalous dimension yet to be established. 

As with all large computations carried out symbolically, it is worth detailing the various 
checks we used in order to be confident that our results are credible. First, for the case of 
the Wilson operators the three loop MS anomalous dimensions are known, [B El O Hj, and 
our anomalous dimensions must therefore agree before extracting any finite parts for lattice 
matching. Moreover, for both Wilson and transversity operators the MS expressions have been 
shown to be independent of the linear covariant gauge fixing parameter. For the transversity case 
we have the checks that the two loop anomalous dimensions must agree with [261 [27] [28"1 [25l [29] 
for the various n we consider. At three loops the only substantial check is that the residues of 
the poles in 1/e 2 and 1/e 3 have to agree with the renormalization group consistency check. In 
other words these are predicted from the one and two loop parts of the anomalous dimensions. 
In addition for the RI' scheme, one can compute the anomalous dimensions either directly from 
the renormalization constants deduced from the Green's function itself, or from the conversion 
functions, Ce>(a, ct), based on the renormalization group. This is defined by 



where the renormalization constants are both expressed in terms of the variables defined in the 
same scheme. Then the RI' anomalous dimension is given by 



(See, for example, [33] ■) Therefore, the expression on the left side will agree with the direct 
renormalization. For all the results presented in the next section we note that they all pass the 
checks discussed here. 






70 (a) - /3(a)— In Cq (a, a) - aj a (a)— In C a (a, a) . 
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6 Results. 



First, we record the explicit values for the anomalous dimensions of the Wilson operators n = 5 
and 6 in RI' for arbitrary a, which are new. The notation is that the numerical subscript denotes 
the moment whilst the superscript, W or T, corresponds to either the Wilson or transversity 
operator respectively. We find 



7T(a) 



and 



91 
15 



+ 



C F a + 



(33525a 2 + 100575a + 1729270) C A 



- 156114C F - 6738807>iV / ] 



C F a 2 



(30172500a 4 + 289359000a 3 



27000 

97200000C 3 a 2 + 1409428125a 2 
- 1004400000C 3 a + 4758071625a - 5142528000C 3 + 52067172425) C\ 
+ (23726250a 3 + 30956400a 2 - 415630950a 

+ 102384000C 3 - 9145680720) C A C F + 60234848007$ iV/ 
- (268200000a 2 - 259200000C 3 a + 1582173000a 

+ 2514240000C 3 + 36792205400) C A T F N f 
+ (107259600a + 3680640000C 3 - 3053173120) C F T F N f 

C F a 3 



+ (1832544000C 3 - 829297168) C F 



48600000 



+ 0(a 4 



(6.1) 



709^ 

C F a + 

105 



(12116475a 2 + 36349425a + 670295290) C A 

C F a 2 



- 57119598C F - 263443880^^] 



+ 



9261000 

(534336547500a 4 + 5228103069000a 3 - 1750329000000C 3 a 2 

+ 25439835416625a 2 - 18086733000000C 3 a + 86004002776125a 
- 92121315552000C 3 + 988839358918775) C\ 
+ (355203339750a 3 + 158333464800a 2 - 9375191062650a 

- 5509035504000C 3 - 172078530172080) C A C F 
- (4749658200000a 2 - 4667544000000C 3 a + 28468726629000a 

+ 49102562880000C 3 + 704961641573000) C A T F N f 
+ (2419404145200a + 73311557760000(3 - 59288998908160) C F T F N f 
+ (31055615136000C 3 - 13674447985168) C\ 

C F a 3 



+ 117065906115200T$A^ 2 



816820200000 



+ 0(a 4 



(6.2) 



where Q n is the Riemann zeta function. As the Landau gauge is of particular interest, we record 
that the previous two expressions when a = are 



a=0 



Q1 C Cr?n 2 
^ F -a + [864635CA-78057CF-336940TFJV/]: 



15 



13500 



+ 



(52067172425 - 5142528000C 3 )Cl 
+ (102384000C 3 - 9145680720)^ C F 
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(2514240000C 3 + 36792205400) C^iV/ 
(1832544000C 3 - 829297168)Cf, 



+ (3680640000C 3 
+ 6023484800T$iV/ 



3053173120)Cjr2>AT / 

C F a 3 , 4- 

+ U(a ) 

48600000 v ; 



(6.3) 



and 



l'e W (a) 



a=0 



709C F 
105 



+ [335147645C A - 28559799C F - 131721940^7^] 
+ 



C F a 2 



4630500 

(988839358918775 - 92121315552000C 3 )Ci 

- (5509035504000C 3 - 172078530172080)C j4 C F 

- (49102562880000C 3 + 704961641573000)C j4 r F iV / 
+ (31055615136000C 3 - 13674447985168)C| 

+ (73311557760000C 3 - 59288998908160)C F T F iV / 

C F a 3 



+ 117065906115200T|iV/ 



+ 0{a q 



(6.4) 



816820200000 

For further comparison between the schemes the MS and RI' expressions for n = 5 transversity 
are, 



ilia) 



92 
15 

+ 



C F a + [189515C A - 41674C F - 798107>iV / ] 

(190836000C 3 + 1975309075) C\ 

- (572508000C 3 + 325464235) C A C F 

- (1192320000C 3 + 511395100) C A T F N f 
+ (381672000C 3 - 254723696) C\ 

+ (1192320000C 3 - 989903260) C F T F N f 
C F a 3 



C F a 2 
6750 



- 83718800T|A/ 



12150000 



+ 0(a 4 



(6.5) 



and 



7 5 



IT i 



92 
— < 
15 



C F a + [(30825a 2 + 92475a + 1740690) C A 
- 166696C F - 676560T F iV / ] 



C F a 2 
27000 



+ 



(l94197500a 4 + 1854279000a 3 - 583200000C 3 a 2 + 8993896875a 2 
- 6026400000C 3 a + 30074295375a - 37353312000C 3 
+ 356401468700) C\ + (91239750a 3 - 209956950a 2 

- 4997987400a + 1076976000C 3 - 60979980560) C A C F 
- (1726200000a 2 - 1555200000C 3 a + 10041363000a 
+ 17858880000C 3 + 253330505600) C A T F N f 



+ 416296832007$ AT/ 

+ 

+ 



(1289803200a + 27164160000C 3 
(IO6868I6OOOC3 - 7132263488) C\ 



22363266560) C F T F N f 
C F a 3 . 4. 

340200000 + ^ ' 



(6.6) 
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For completeness we record the next MS anomalous dimensions in the sequence are, [20J, 

C F a 2 



I 1 , (a) 



and 



77 T («) 



34 „ 

—C F a + 

5 



[204770C A - 42129C F - 888WT F N f ] 



6750 



+ 



(707616000C 3 + 7527909825) C\ 

- (2122848000C 3 + 1373507730) C A C F 

- (4626720000C3 + 1841332000) C A T F N } 
+ (1415232000C 3 - 684744816) C| 

+ (4626720000C 3 - 3910683210) C F T F N } 

C F a 3 



3209758007$ AT/ 



42525000 



+ 0(a q 



(6.7) 



258 „ 

— C F a + 
35 



[75266555CA - 15484767Cf - 331498307>iV/] 



C F a 2 
2315250 



+ 



(3517994592000C 3 + 38365845513450) C\ 

- (10553983776000C 3 + 5978407701105) C A C F 

- (24084527040000C 3 + 9039144860900) C A T F N f 
(7035989184000(3 - 4192441946262) C| 
(24084527040000C 3 - 20698675427220) C F T F N f 

C F a 3 



+ 
+ 



16513111916007$ AT/ 



+ 0(a 4 



(6.8) 



204205050000 

The complete set of three loop transversity anomalous dimensions in MS and RI' are given in 

[Tsunami up. 



7 Conclusions. 



We conclude with a few brief remarks. First, the three loop anomalous dimensions are available 
for the transversity operator for each moment up to n = 8 in the MS scheme and n = 7 for 
the lattice motivated RI' scheme. The former in particular will provide important independent 
checks for future explicit arbitrary moment evaluations of the three loop anomalous dimensions. 
A by-product of the overall project, [THJ [191 I2Q], has been the provision of the finite parts 
of a Green's function which are necessary for lattice measurements of matrix elements. The 
three loop perturbative information is essential to obtaining more precise numerical estimates. 
In addition we have given the (new) RI' anomalous dimensions for the n = 5 and 6 Wilson 
operators at three loops. Whilst it is in principle possible to continue with the computation of the 
transversity higher moments to n = 9 and beyond, the present method has become too tedious. 
This is primarily due to the increase in the number of free Lorentz indices on the operator 
which was originally required for the lattice comparison. Moreover, the actual computation 
time as indicated for n = 8 in the Feynman gauge has already become unacceptably long. An 
explicit arbitrary n calculation exploiting the algorithm of [U [21 S] would achieve all moment 
information during one run, at possibly a computation time which is not too dissimilar from 
that for one high moment. 
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